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ABSTRACT 


'*  A  satellite  has  been  observed  through  an  optical  telescope 
with  a  CCD  camera,  the  electrical  signal  stored,  and  later 
analyzed  with  the  aid  of  a  digital  computer.  The  unresolved 
satellite  image  is  contained  in  3  to  7  picture  elements  in  fifty 
contiguous  data  frames.  This  allows  a  better  understanding  of 
the  physical  characteristics  of  the  satellite  image  as  viewed 
through  the  atmosphere  by  an  electro-optical  detector. 

This  paper  describes  the  data  processing  programs  used  to 
analyze  the  data,  the  analysis  results,  and  implications  for 
signal  processing  design. 
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I .  INTRODUCTION 

A  CCD  camera  has  the  potential  of  improving  the  performance 
of  satellite  detection  systems  by  means  of  higher  quantum 
efficiency,  smaller  pixel  size,  and  more  pixels  on  a  focal  plane 
than  EBSICON  cameras.  An  efficient  satellite  detection  processor 
should  be  designed  to  exploit  the  signal  characteristics  of  the 
camera  it  uses  as  an  electro-optical  detector.  Selected  portions 
of  a  digital  satellite  detection  processor  (IMC  Signal  Processor) 
and  CCD  camera  are  being  developed  for  the  GEODSS  advanced  tech¬ 
nology  program.  The  output  of  a  CCD  camera,  when  observing  a 
satellite,  has  been  digitized,  and  then  analyzed  with  the  aid  of 
a  computer.  This  allowed  a  detailed  study  of  the  electrical  signal 
that  the  processor  would  have  to  use  prior  to  designing  the 
processor . 

Satellite  ATS-6  was  observed  with  the  31"  main  GEODSS  ETS 
telescope  at  the  White  Sands  Missile  Range.  The  CCD  camera  was 
at  the  prime  focus  and  was  operated  in  the  stare  mode  with  an  0.1 
second  integration  time.  The  ETS  photometer  measured  the 
satellite's  brightness  as  10.5  Mv.  The  camera  video  was 
digitized  to  12  bit  resolution,  sent  to  the  MODCOMP  computer 
at  the  ETS  and  stored  on  magnetic  tape.  This  tape  was  brought 
back  to  Lexington  and  placed  on  the  MODCOMP  computer  here.  The 
tape  was  read  by  the  MODCOMP  and  the  data  transmitted  to  the 
HP  9845  graphics  computer  and  stored  there.  The  HP  computer 
was  used  for  all  the  data  analysis. 
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The  satellite  observation  is  contained  in  fifty  contiguous 
frames.  Each  frame  contains  an  array  of  32  x  32  picture  elements 
(pixels).  The  integration  time  for  each  frame  was  0.1  second. 
Each  pixel  has  the  equivalent  field-of-view  of  2” 7  x  2" 7. 


II.  DATA  PROCESSING 

The  frames  of  data  were  read  from  the  digital  tape  by  the 

MODCOMP  computer  and  transmitted,  on  an  RS-232C  interface1,  to 

.  2 
the  HP  computer.  Using  this  interface  the  HP  computer  was 

operated  to  look  like  a  standard  terminal  to  the  MODCOMP.  The 

major  task  of  this  interface  was  to  write  an  HP  terminal 

emulator  program.  This  program  has  to  read  the  data  from  the 

interface,  store  the  data,  and  respond  to  the  MODCOMP  with 

standard  terminal  handshake  commands.  The  program,  "MOD-HP", 

listed  in  Appendix  A  did  this  at  9600  baud. 

To  operate  at  this  high  a  rate,  each  data  word  was  read 

as  a  character  string.  The  character  strings  were  stored  on 

HP  tape  cassettes  as  the  data  word  was  being  read.  The  1024 

character  strings  from  each  frame  were  separated  by  markers. 

Once  all  fifty  frames  of  the  character  strings  were  stored  on 

tape,  the  HP  computer  was  disconnected  from  the  MODCOMP,  and 

was  operated  independently  from  then  on. 

Two  tape  cassettes  were  required  to  store  all  the  character 

strings.  All  the  character  strings  were  then  read  off  the  tapes 

and  stored  on  one  floppy  disc  using  the  "COPY"  program  listed 

in  Appendix  A.  The  character  strings  cannot  be  used  in 

programs,  but  must  be  converted  to  the  numbers  they  represent. 

The  program  "READ",  Appendix  A,  read  the  character  strings  off 

the  disc,  converted  the  character  strings  to  numbers  and  stored 

l"!  Hewlett  Packard,  "Serial  I/O  Interface-98036A, "  1977. 

2.  Hewlett  Packard,  "I/O  ROM  Programming  System  45B,"  1979. 


these  on  tape.  The  "COPY"  program  was  again  used  to  transfer 
data  from  tape  to  disc. 

To  insure  that  the  data  stored  on  the  disc  was  accurate, 
a  print-out  of  the  data  stored  on  the  MODCOMP  tapes  was  compared 
to  the  data  stored  on  the  HP  disc.  The  data  was  identical  in 
frame- to- frame  and  entry-to-entry  comparisons.  The  first  data 
entry  in  each  frame  of  both  the  MODCOMP  and  HP  data  was  the 
same,  but  was  not  a  valid  data  value.  To  avoid  problems  in 
data  processing  the  value  of  the  first  data  entry,  in  each  frame, 
was  set  to  the  mean  of/  the  data  of  each  frame  and  is  then  called 
a  valid  data  frame. 

Now  the  data  is  ready  to  be  analyzed.  The  first  data 
analysis  program  is  "MAXAMP"  listed  in  Appendix  B.  "MAXAMP”  finds 
the  maximum  amplitude  in  each  frame  and  prints  the  frame  number, 
amplitude,  and  array  position.  A  sample  of  "MAXAMP"  output  is 
shown  in  Figure  1.  The  frame  number,  e.g.  T11V01:F  stands  for 
CCD- DIGITIZER  test  #11  (Til) ,  valid  data  frame  (V) ,  frame  number 
(01) ,  and  stored  on  floppy  disc  (F) .  The  pixel  array  position 
is  numbered  from  1  to  1024  in  each  data  frame. 

The  next  data  analysis  program  listed  in  Appendix  B  is 
"FREDIS "  which  plots  the  frequency  distribution  of  amplitudes 
for  each  frame.  Figure  2  contains  plots  from  two  of  the  fifty 
frames.  The  amplitude  scale  adjusts  for  each  frame  depending  on 
the  maximum  amplitude.  The  occurrence  scale  also  adjusts 


4 


fill 


FRAME 

T  1  1  V  1  0 :  F 

MAXIMUM 

i  1 1  v  1 1 :  f 

MAX  I  MUM 

T 11 V 1 2 : F 

MAXIMUM 

T  1  i  V  1  3  :  F 

MAXIMUM 

T 1 IV 1 4  :  F 

MAXIMUM 

T  1 1  V  1  5 :  F 

MAXIMUM 

T  i  1  V  1  6 :  F 

MAX  I  MUM 

^  1  1  V  1  7  :  F 

MAXIMUM 

T  I  1  V  1  8 :  F 

MAXIMUM 

T i 1 V 19 : F 

MAXIMUM 

T 1 1728: F 

MAXIMUM 

11 1721 : F 

MAXIMUM 

T 1 1722: F 

MAXIMUM 

r i 1723: f 

MAXIMUM 

r  1 1724:  f 

MAXIMUM 

T i 1725: F 

MAXIMUM 

Fig.  1. 


AMPLITUDE 

185 

AMPLITUDE 

182 

AMPLITUDE 

1 6  3 

AMPLITUDE 

159 

AMPLITUDE 

224 

AMPLITUDE 

176 

AMPLITUDE 

208 

AMPLITUDE 

173 

AMPLITUDE 

147 

AMPLITUDE 

1  80 

AMPLITUDE 

172 

AMPLITUDE 

143 

AMPLITUDE 

144 

AMPLITUDE 

143 

AMPLITUDE 

175 

AMPLITUDE 

134 

"MAXAMP"  sample  output. 


5 


POS  I  T I  ON 


i_i9  I 


depending  on  the  number  of  pixels,  excluding  the  zero  amplitude, 
for  which  there  are  from  350-400  pixels  in  each  frame. 

The  program  "PEAKS"  examines  a  frame  of  data  and  prints  out 
the  location  and  amplitude  of  any  pixel  whose  amplitude  is  above 
a  selected  level.  Figure  3  is  a  sample  of  "PEAKS"  output.  A 
program  listing  of  "PEAKS"  is  in  Appendix  B.  The  frame  number 
is  given  and  the  selected  level  (e.g.  19).  The  row,  column, 
array  position  (as  in  "MAXAMP"  output) ,  and  corresponding 
amplitude  are  listed. 

Program  "DATAST",  Appendix  B,  prints  out  the  amplitude 
values  of  any  frame  by  row  and  column.  An  asterisk  is  placed 
next  to  amplitudes  above  a  selected  amplitude.  The  output  is 
displayed  on  the  CRT  or  printed  on  the  line  printer.  A  sample  of 
output  from  frame  #11  is  shown  in  Figure  4.  The  top  line  of 
each  pair  is  array  elements  1  to  16  and  the  second  line  is 
array  elements  17  to  32.  In  this  frame  the  satellite  pixels 
have  amplitudes  of  182,  112,  32,  and  40.  An  isolated  pixel  of 
noise  with  amplitude  32  is  one  array  line  above  this. 

The  "FRAME"  program  (Appendix  B)  allows  a  13  x  11  pixel 
array  around  the  satellite  image  of  any  frame  to  be  displayed. 
After  the  array  of  pixels  is  displayed,  the  amplitudes  are 
plotted  first  in  column  and  then  in  row.  Figure  5  is  a  print¬ 
out  of  the  array  of  pixels  around  the  satellite  image  from 
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frame  #12.  Figure  6  is  a  plot  of  the  amplitude  in  the  sixth 
column  of  Figure  5.  The  dashed  line  is  a  level  at  an  amplitude 
of  19.  Figure  7  is  a  plot  of  amplitudes  in  the  fifth  row  of 
Figure  5.  The  dashed  line  is  also  at  an  amplitude  of  19. 

Appendix  C  contains  the  pixels  (5  x  6)  surrounding  the 
satellite  image  for  frames  1  through  12.  These  were  obtained 
using  a  modified  version  of  program  "FRAME" . 
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III.  ANALYSIS  OF  DATA 


Analysis  of  the  output  of  the  "MAXAMP"  .program  indicates 
that  the  highest  amplitude  of  the  satellite  image  moves  between 
two  adjacent  pixels.  Over  the  fifty  frames,  the  highest 
amplitude  switches  15  times  between  the  two  pixels ,  spends  an 
equal  amount  of  time  in  each  and  there  is  no  periodicity  to  the 
switches.  The  range  of  highest  amplitudes  is  127  to  224  with 
a  mean  of  165.  The  mean  of  highest  amplitude  in  one  pixel  is 
162,  while  the  other  pixel's  highest  amplitude  mean  is  169. 

Examining  the  output  of  "FREDIS"  shows  that  in  all  frames 
most  of  the  pixels  have  amplitudes  from  0  to  17.  The  mean 
of  the  amplitudes  from  each  of  fifty  frames  ranges  from  4.4  to 
4.9.  This  is  strongly  due  to  the  large  number  of  pixels  (350  to 
400)  in  each  frame  with  amplitudes  of  zero.  The  distributions 
of  all  fifty  frames  are  similar.  All  frames  have  less  than  20 
pixels  above  an  amplitude  of  19. 

The  output  of  "PEAKS"  shows  that  in  every  frame  the  pixels 
with  amplitudes  greater  than  28  are  either  those  of  the  satellite 
image  or  of  isolated  pixels.  There  are  532  pixels  that  have 
amplitudes  greater  than  19  in  the  fifty  frames.  Most  of  these 
pixels  are  either  those  of  the  satellite  image  or  isolated 
pixels.  However,  12  of  the  pixels  in  10  of  the  fifty  frames 
have  an  adjacent  pixel  greater  than  19.  In  half  of  these  cases 
the  adjacent  pixels  are  vertically  related  and  in  half  horizontally. 


The  pixels  with  amplitudes  greater  than  39  are  all  those  of 
the  satellite  image. 

The  output  of  "DATAST"  shows  that  the  noise  is  uncorrelated 
from  pixel  to  pixel  and  from  line  to  line.  Eight-four  percent 
of  the  pixels  with  an  amplitude  greater  than  8  have  no  adjacent 
pixel  with  amplitude  greater  than  8.  Ninety-three  percent  of 
the  pixels  with  amplitude  greater  than  12  have  no  adjacent  pixel 
with  amplitude  greater  than  12. 

The  output  of  "FRAME"  gives  information  primarily  about  the 
pixels  associated  with  the  satellite  image.  Table  1  is  a  summary 
of  the  distribution  of  satellite  image  pixels  by  frames.  The 
distribution  of  energy  from  the  satellite  image  is  very  dynamic. 
The  distribution  changes  significantly  in  periods  of  0.1  sec  to 
0.3  sec.  In  61.5%  of  the  frames  the  energy  distribution  doesn't 
fit  a  circular  Gaussian  function. 
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TABLE  1 


NUMBER  OF  SATELLITE  PIXELS  FOR  FIFTY  FRAMES 


Amplitude  Higher  Than 
24 


Satellite  Pixels 

3 

4 

5 

6 
7 


Number  of  Frames 
10 
25 
12 
2 
1 


29 


2 

3 

4 

5 


1 

19 

28 

2 


32 


2 

3 

4 


17 

31 
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34 
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39 


2  29 

3  21 
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44 
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99 


1  16 

2  34 
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IV.  SIGNAL  PROCESSING  CONSIDERATIONS 

After  reviewing  the  analysis  results,  it  is  apparent  that 
the  satellite  image  is  not  stable  over  time  periods  as  short 
as  0.2  sec.  The  instability  of  the  image  could  result  from  three 
primary  sources:  satellite  motion,  telescope  jitter,  and  atmospheric 
seeing.  During  the  time  of  data  collection,  the  ATS-6  satellite 
was  a  three-axis  stabilized  geosynchronous  satellite.  Its 
reflected  light  vs.  time  had  been  measured  by  several  observers 
at  the  ETS  and  was  found  to  be  steady  throughout  the  night,  as 
would  be  expected  from  this  type  of  satellite.  Therefore  satellite 
motion  didn't  account  for  the  image  instability  observed  in  the 
data.  The  fifty  frames  of  data  were  carefully  scrutinized  to 
detect  any  indication  of  telescope  jitter.  Telescope  jitter 
would  be  detected  in  the  data,  if  the  pixels  which  contain  the 
satellite  image  shifted  together  in  the  same  direction.  There  is 
no  observable  correlation  between  the  maximum  amplitude  pixel  shift 
direction  and  the  other  satellite  image  pixels'  shift  direction. 
Calculating  the  centroid  of  the  satellite  pixels  shows  that  the 
centroid  movement  is  in  random  directions  rather  than  in  one  or 
two  orthogonal  directions  that  would  characterize  telescope  jitter. 
Atmospheric  seeing  is  therefore  the  leading  candidate  as  the 
source  of  the  image  instability.  Although  few  measurements  of 
atmospheric  seeing  with  time  constants  as  short  as  0.1-0. 2  sec, 
and  with  two-dimensional  spatial  resolution  as  fine  as  277  have 
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been  reported,  one  reference  does  report  the  characteristic 
image  fluctuation  seen  in  this  data  as  atmospheric  seeing.  Many 
other  reports  on  atmospheric  seeing  describe  image  fluctuation 
similar  to  that  observed  in  this  data.  Atmospheric  seeing  will 
affect  stellar  images  in  the  same  fashion.  A  signal  processing 
system  that  removes  stars  by  comparing  a  star  field  at  one  time 
with  the  same  star  field  at  a  time  0.2  sec  later  or  longer,  must 
incorporate  a  technique  to  reduce  false  alarms  due  to  the  star 
image  instability.  These  false  alarms  are  located  adjacent  to 
star  images,  and  a  guard  band  around  each  star  that  inhibits 
detections  would  remove  these  false  alarms. 

Although  the  resolution  of  the  telescope  is  1'.'0  or  smaller 
and  the  pixel  resolution  is  2'. '74,  the  satellite  image  is  in 
3  to  7  pixels.  Most  of  the  energy  is  in  2  to  3  pixels,  however 
in  74%  of  the  frames  the  energy  is  in  4  to  5  pixels.  In  more 
than  half  the  frames  the  distribution  of  satellite  energy 
doesn't  fit  into  a  circular  function.  There  is  little  correla¬ 
tion  from  frame-to-f rame  of  the  spatial  distribution  of  satellite 
energy.  This  requires  that  predictions  of  satellite  SNR  and 
system  detection  performance  be  modified  to  incorporate  these 
physical  realities. 


3.  MI RESELL,  A.H.  et  al,  "The  Scintillation  of  Starlight",  JOSA, 
41,  689-695,  1951. 
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The  star  and  satellite  images  are  in  several  adjacent  pixels 
(at  a  2'.'74  pixel  size),  while  the  noise  is  only  one  pixel  long. 
This  is  an  excellent  discriminant  between  the  satellite  and  high 
amplitude  noise  that  competes  with  the  satellite  signal.  If 
a  threshold  level  of  28  was  used  with  the  pixels  in  the  fifty 
frames,  the  threshold  crossings  would  be  due  to  the  satellite 
(adjacent  threshold  crossings)  or  to  the  noise  (isolated  thres¬ 
hold  crossings).  The  isolated  threshold  crossings  could  then  be 
removed  and  only  the  satellite  threshold  crossings  would  remain. 
As  the  threshold  is  lowered,  noise  threshold  crossings  will  be 
more  frequent  and  false  alarms  due  to  noise  will  start  to  re¬ 
appear.  At  a  threshold  of  28,  100%  of  the  noise  false  alarms 
are  removed;  at  a  threshold  of  12,  93%;  and  at  a  threshold  of 
8,  84%.  The  remaining  false  alarms  due  to  noise  can  be  removed 
by  additional  processing  (streak,  slope,  velocity,  etc.). 
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V. 


CONCLUSIONS 


The  analysis  of  a  satellite  image  as  seen  by  a  CCD  camera 
mounted  on  the  31"  prime  focus  at  the  ETS  has  been  accomplished. 
The  complexity  of  the  image  as  a  function  of  spatial  distribution 
and  time  can  be  clearly  seen.  The  seeing  effects  were  worse  than 
what  had  previously  been  expected  at  the  ETS.  The  utility  of 
having  a  digitized,  stored  signal  from  a  sensor  before  designing 
a  signal  processor  that  uses  that  particular  sensor,  is  obvious. 

The  pixel  size  of  a  camera  is  important  to  the  detection 
system;  as  too  small  a  pixel  will  cause  the  satellite  energv  to 
be  spread  over  many  pixels  and  thus  reduce  the  SNR,  or  if  the 
pixel  is  too  big  the  night  sky  will  contribute  additional  noise 
and  also  lower  the  SNR.  The  image  size  depends  on  the  atmospheric 
seeing,  therefore,  the  selection  of  pixel  size  depends  on  the 
seeing  expected  at  a  site.  The  ideal  way  to  measure  seeing  at 
the  ETS  would  be  to  use  an  automated  version  of  what  was  done 
in  this  report.  The  microprocessor  based  PVDS  system  would  be 
ideal  for  this  purpose  and  could  provide  a  short  exposure 
measure  of  seeing.  Using  a  CCD  camera  mounted  on  the  cassegrain 
position  of  the  31"  ETS  telescope  would  provide  even  more  detail 
of  the  image  intensity  distribution.  Binary  stars  could  be 
viewed  to  eliminate  some  observing  uncertainties  and  data  could 
be  collected  on  a  regular  basis. 
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APPENDIX  A 


1  0 

DIM  D*'.  500  .>[362] 

!  " MOD-HP “ 

20 

l REATE  “TtlFe  I'M,  10000 

30 

ASSIGN  #1  TO  " T i 1 F 0 1 “ 

40 

RESET  10 

50 

WRITE  10  10,5; 1 

60 

WAIT  WRITE  10,4564 

70 

WAIT  WRITE  10, 4; 203 

SO 

WAIT  WRITE  10,4; 32+7 

90 

READ  10  10,4;R4 

100 

WRITE  10  10, 5; 0 

1  1  u 

WAIT  READ  10, 4; A 

! READY 

1  20 

WRITE  10  10. 7; 0 

! TRIGGER 

13  O  i 

14o  *  SET  UP  THE  98036*  TO  INTERRUPT  ON  INPUT  DATA 

150 

1*0  WRITE  10  10,5; 128+4 

i ro  p  =  i oflag • i 0  > 

1  80  8  =  1  OSThTUSC  1  0  1 

190  PRINTER  IS  1* 

100  PRINT  "I/O  FLAG  1  "  ,  F  ,  "  I  -  0  STATUS  1  " , S 
210  READ  10  10,5;R5 

220  READ  1 0  10,*; P* 

230  STATUS  10; V 

240  PRINT  " R4  0  "  ,  R4 , " R5  146  ",R5,"R6  248  *' ,R€, "V  402  "  ,  V 

250  1=0  ‘INDEX  FOP  WORDS  COMING  OVER 

260  BEEP 
270  BEEP 

2 8 0  T 0 P E N  10,2  G 0 S U B  Get i n p u t 

290  Loop:  !  WRIT  FOR  DATA  FROM  HODCOMF' 

300  ON  KBD  1  GOTO  310 
310  K t  =  K  B  D  t 

320  Ex*POSCKT,  CHR*  t' 255>««CHR$< .21  >  > 

3  30  IF  E  s<  >0  THEN  440 

340  GOTO  Loop 

350  Get  input:  !  GET  INPUT  FROM  INTERFACE 

360  1=1+1 

370  IF  I  5O0  THEN  Loop 
380  D*< I > =TBUF $ 

390  IF  P  0  S  D  %  '  I  >  ,  C  H  R  $  <  5  >  >  THEN  OUTPUT  10;  CHR*  (.  6  > 

400  PRINTER  IS  16 

410  PRINT  I 

420  WHIT  900 

430  RETURN 

440  FOR  1=1  TO  500 

450  PRINT  #1 ; I , b$<  I  > 

460  NEXT  I 
470  END 


"MOD-HP* 
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FOR  H=10  TO  50  !  "feEAD" 

D  I  M  DS  C  2O0  >  C  320  ]  ,  F‘  '■  32  1  >  ,  0  <  200  ,  F  C  1 024  > 

X  S  =  V  H  L  S  1  W  > 

W1S="T1 1F“£XS$ M:F8M 
W2*="T1  1  D"  &XS&  M  :  T 1 5  *’ 

ASSIGN  #1  TO  W1S 
ON  END  #1  GOTO  660 
St  —  0 
D  =  0 

FOP  1-1  TO  200  !  READ  STRING  FROM  TAPE  C  I' 

READ  #1 ;  I  ,D*C I > 

DS  <  I  >  =  DS  <1  >  C  1 , 3 1 9  3 &CHRS  •-  1 3  > 

PRINTER  IS  16 

FOR  J=1  TO  200  !  BREAK  STRING  INTO  VARIABLE 

IF  J ; 1  THEN  290 
F  C  1  =POS  (  DS  <  I  >  ,  CHRS  Cl  3  >  > 

Lf =POS < DSC  I > 1 1 , P Cl > 3 , CHRS C 1 0 > > 

L  =  LEN'  DSC  I  > C 1 , PC  1  > 3 > 

IF  Lf < >0  THEN  240 

IF  L< 2  THEN  280 

OC  1  >  =  VAL C  DSC  I  >C  1  ,  PC  1  >  ]  :> 

D=1 

GOTO  290 

IF  L<3  THEN  280 

0  <  1  >  s  V  A  L  C  D  S  C  I  >  C  2 ,  P  C  1.0) 

D  =  1 

GOTO  290 
OC  1  >=-l 

ji«j+i 

P  C J 1 >  =POS CDSC I )[PCJ)  +  1 3, C HR SC  1 3  > > 

IF  PCJ1>=0  THEN  480 
PC  J1  >=PC  J1  >+PC  J> 

L=LENCDSC I >CPC J)+l , PC J1 > ] > 

IF  L< 4  THEN  460 
Lf  =POS CDSC I >, CHRS C 1 0 >  > 

IF  Lf < >0  THEN  390 
G  =  1 

GOTO  430 

lj  =  2 

DS C  I  >  C P C J >+G, PC  J 1 >  3-TRI MS C  DS  U  > C P < J  >  +G , P C J 1 >  ]  > 
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4  1  Ci 

420 
4  30 
440 
450 
460 

4  ro 

480 

480 

500 

510 

520 

530 

540 

550 

5  6  0 
570 
530 
580 
600 

6  1  u 
620 
6  3  O 
6  4  0 
6  5  0 
660 
670 
6  3  0 
6  80 
700 
710 
720 
730 
740 
750 
760 
770 
780 
790 
300 


L=LEN<  D*<  I >[P<J>+G,P<Jl>3> 

IF  L  6  THEN  460 

0 < J 1 >  =  V  A  L (D*<I>CP<J>+G,P(J1>1) 

B  =  D  +  1 

GOTO  470 

D< J1  >*-l 

GOTO  5O0 

S  =  D 

GOTO  510 
NEXT  J 


35  =  0 

C  =•  =  S  t  +  1  -  S  s 
C  f  =  3 1 +s 
D  =  0 

IF  < D < 1 >  =  - 1 >  ANT  <S*0>  THEN  670 
F 0 R  K  =  C 5  TO  C f 
L  =  K  -  S  t  +  $  s 

IF  0 v  L  >  = - 1  THEN  600 

GOTO  640 

L  =  L+  1 

C  f  =  C  f  - 1 

Ss  =  Ss  + 1 

GOTO  580 

F  <.  K  >  =  0  <  L  > 

NEXT  K 
St =St +s 
NEXT  I 

FOR  M=1  TO  32 

FOR  N=M  TO  102  3  +  M  STEP  32 

IMAGE  #,4B 

PRINT  USING  700;FCN> 

NEXT  N 

PRINT  L I N < 1 > 

NEXT  M 


ON  END  #2  GOTO  790 
CREATE  W2$, 1,9000 
ASSIGN  #2  TO  W2* 


PRINT 

#2;f<*> 

NEXT 

14 

END 

10 

FOR  1=10  TO  50 

20 

I  $  =  V  A  L  ♦  <  I  > 

30 

T  $  = " T 1  ID” 

40 

Pit=H : T15" 

50 

P2$=" : F" 

60 

F 1  $*T$& 1  f&P  1  ♦ 

70 

F2*  =  T*U*&P2  % 

80 

PRINT  F 1  $ , F2$ 

90 

COPY  Fll  TO  F2* 

100 

NEXT  I 

110 

END 

"  READ” 


(Cont.) 


NEXT  J 


PACK  VARIABLES  FOR  OUTPUT  K 


NEXT  K 
NEXT  1 

OUTPUT  DATA  <M> 


NEXT  M 


!  "COPY" 

!  FRAME  NUMBER 

!  TAPE  DRIVE 
!  DISC 


"COPY* 
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10  PRINT  "FRAME","  ","  "."POSITION"  !  "MAX AMP" 

20  FOR  1=10  TO  50 

30  DIM  F <1 024) 

40  I*=VAL*< I > 

50  V*="TllV"«iI*S,":F" 

60  ASSIGN  #1  TO  V* 

70  ON  END  #1  GOTO  100 

80  READ  # 1 ; F  <  *  > 

90  Max =8 

100  FOR  M= 1  TO  32  !  OUTPUT  DATA  <M> 

110  FOR  N=M  TO  1023+M  STEP  32 

120  IF  F ( N ) >Max  THEN  140 

130  GOTO  160 

140  M  ax  -  F  (  N  ) 

150  Na:;=N 

160  NEXT  N 

170  NEXT  M  !  NEXT  M 

ISO  PRINT  V*, "MAXIMUM  AMPL I TUDE " , Max , Nax 
190  PRINT  LIN<1 > 

200  NEXT  I 

210  END 


"MAXAMP" 


10  PRINTER  IS  16  !  H FRED  1 9" 

£0  PRINT  PAGE;"  This  program  allows  the  user  to  enter  a  frame  of  data  and  ha 

we  a" 

30  PRINT  "  frequency  distribution  computed  and  drawn." 

40  PRINT  LINO);"  After  the  distribution,  a  beep  will  sound.  The  chart  will 

50  PRINT  "remain  on  the  screen  and  the  program  will  halt  until  you  press  CONT 

60  PRINT  "This  will  give  you  time  to  digest  the  information  on  the  screen." 

70  PRINT  LINO);"  If  you  want  the  output  of  this  program  dumped  from  the  CRT 

8 0  P  R I N T  " t  o  t h e  i n t  e  r n  a 1  the r  m a 1  printer,  then  an  s  w  e  r  Y  to  the  q u e  s  t i o n  b e 1 o w 

90  PRINT  "If  you  don't  want  hardcopy  answer  N.  Then  press  CONT." 

100  INPUT  "Do  you  want  hardcopy  <Y^N)?",At 

110  IF  U  P  C  $  <  A  $  >  =  "  Y  "  THEN  Yhard 

120  IF  UPC$< A$ )=M N"  THEN  Nhard 

130  BEEP 

140  GOTO  100 

150  Yhard S  Hardcopy* 1 

160  GOTO  130 

170  Nhard!  Hardcopy=0 

130  N  =  224 

190  CALL  Dr i wer< N, Hardcopy > 

20O  END 

210  SUB  Dr i wer<Nf Hardcopy > 

220  RANDOMIZE 

230  DIM  A < N > , L$ < N ) , Un i t s* C30 3 , V< 1 024 ) 

"FREDIS" 
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24©  INPUT  "Which  frame  01  to  5©?",Fr* 

250  V$  =  "Til  V  "  &Fr  ** :  F  " 

26©  ASSIGN  #1  TO  V$ 

2?0  REAI»  #l;V<*> 

200  Max- V  < 1 > 

290  FOR  J= 1  TO  1024 
300  Max=NAX < Max , V < J >  > 

310  NEXT  J 
-320  MAT  A = Z  E  R 

3  30  H i g h  =  O 

340  Lh i gh=4©96 
350  FOR  K*  1  TO  1024 

360  I F  ( V < K  >  >H i gh >  AND  < V i K >  <  Lh i gh >  THEN  380 
370  GOTO  390 
3 SO  H l g  h  =  V \ K > 

390  NEXT  K 

40O  FOR  Jl  =  l  TO  1024 

410  IF  V < J 1 > s H i g h  THEN  430 

420  GOTO  440 

4 30  A  '•  H  i  gh  >  =  A  v  H  i  gh  >  +  1 

440  NEXT  J1 

450  IF  High*©  THEN  490 

4b0  L h l g  h = H l g h 

470  High*© 

4  3  O  G  Ci  T  0  3  5  O 
490  PRINTER  IS  © 

5O0  PRINT 

510  PRINT  V* 

520  PRINT  "MAX  AMP" 

530  PRINT  Max 
540  PRINTER  IS  16 

550  FOR  I=20  TO  Max  STEP  2©  !  ENTER  THE  TITLES 

560  L$a>=VAL*a> 

57©  NEXT  I 

5 8 u  U n i t  sl= " N u m b e r  o f  P i x els" 

590  T i t 1 e$* " Amp 1 i  t  ude s “ 

60O  N  =  M  ax 

6 1  ©  CALL  B ar < A  <  *  > , L  $  < *  > , N , U n i t  s  $ , T i 1 1 e  $ > 

620  BEEP 
6 30  PAUSE 

64©  IF  Hardcopy* 1  THEN  BUMP  GRAPHICS 
650  EXIT  GRAPHICS 
660  SUBEXIT 

670  SUB  Bar  ,  L$<*>  ,  N,  Uni  t  »$, Ti 1 1  e$> 

63©  PLOTTER  IS  "GRAPHICS" 

690  GCLEAR 
70O  GRAPHICS 
710  PEN  1 
7  2  0  S  U  Hi  =  0 

"FREDIS"  (Cont.) 
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7  30  Ma>:  =  fid> 

740  FOR  1=1  TO  N 
750  S  L4  m = S  u  m  +  ft  <  I  > 

760  Max  =  M RX  < M ax , R  < I > > 

770  NEXT  I 
700  Ly=LGT<Max) 

790  SETGU 

30O  LOCRTE  15,120,40,95  !  Set  screen  boundaries  and  scale 

810  CLIP  0, 123,40, 100 
820  MOVE  61.5,97 
830  LORG  5 

840  LRBEL  USING  850;Tit1e$ 

850  IMRGE  # , K 

360  SCRLE  0, N+. 5, 0, Max+Max/10 

870  Xtic=0  !  No  tic  marks  on  horizontal  axis 

880  Testyt i c =FRRCT < Ly > + < Ly < 0 > 

3  9  0  V  t  i  c  =  1 0  dNT<Ly>-l)*d  +  1.5*<<Testytic>.39?94  >  AND  (Tes  t  y  t  i  c  <. .  6 9  8  9  7  >  >  +  4  *  <  <  T  € 
st yt i c  >= . 69897)  RND  < Testyt i  c<  =  . 87506) ) +6. 5*< Test yt i c >. 87506) > 

900  !  Compute  the  vertical  tic  marks 

910  ChLL  LaxesCXt i c , Yt i c , 0, 0, 1 , 1 , 2 , 0 , N  + . 5 , 0 , Max+Max/10> 

920  LORG  4 
930  FIXED  1 
940  FOR  1=0  TO  N 
950  MOVE  I -.5,0 

960  I  PLOT  0 , R< I ) , - 1  !  Draw  the  bar  for  the  ith  data  item 

970  I PLOT  1,0 
980  IPLOT  0 ,  -  fi  d  ) 

990  MOVE  I,  fid) 

1W00  NEXT  I 
1010  DEG 
1020  IMfiGE  K 

1O30  LDIR  90  !  Rotate  the  label  90  degrees 

1040  SETGU 
1050  LORG  5 
1060  MOVE  5,70 

1070  LRBEL  USING  1020;  Unit  s$  !  Write  the  d  i  s  p  1  ay  u  n  its  t<  e  i  n  g 

1080  !  used  along  the  vertical  axis 

1090  CLIP  0,123,0,100 

1100  SETUU 

1110  FUR  1=1  TO  N 


"FREDIS"  (Cont.) 
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1120  MOVE  I, -Max  20 
1130  LORG  8 
1 140  PEN  1 

1158  LDIR  90  !  Rotate  the  labels  90  degrees 

1160  LABEL  USING  1020;L*<I> 

1170  NEXT  I 
1180  BEEP 
1190  SUBEXIT 

1 200  SUB  L  axes  <  Xt  i  c  ,  Yt  ic,  Xorg,  Yorg,  Xma  j  ,  Yma.j  ,  M  i  nt  i  c  s  i  zc  f  Xm  i  n,  Xrnax ,  Ym  i  n ,  Ymax  > 
1218  IF  <  X  m  i  n  >  =  X  m  ax  >  OR  <  Y  m  i  n  > = Ymax  >  THEN  S  U  B  E  X I  T 
1 220  Xf udge= . 02*  < Xmax-Xm i n ) 

1 230  Yf udge= . 82* < Ymax- Yin i n > 

1 240  AXES  Xt i c , Yt i c , Xorg, Yorg , Xmaj , Ymaj , M i m i c s i ze 
1258  DEG 

1260  IF  NOT  Xt i c  THEN  Labely 

1278  Labe  1 xt IF  SGN < Xt i c ) *- 1  THEN  Parx 

1280  LDIR  90 

1290  LORG  8 

1300  GOTO  1330 

1310  Parx:  LDIR  8 

1328  LORG  6 

1330  FOR  I=Xorg  TO  Xmax  STEP  ABS<Xtic) 

1340  MOVE  I , Yorg-Yfudge 
1358  LABEL  USING  1360; I 
1360  IMAGE  #,K 
1370  NEXT  I 

1388  Labe  1 y :  IF  NOT  Yt i c  THEN  SUBEXIT 

1398  IF  SGN<Yt i c >*-l  THEN  Pary 

1400  LDIR  0 

1418  LORG  8 

1420  GOTO  1458 

1430  Pary:  LDIR  -90 

1448  LORG  6 

1450  FOR  I = Yorg  TO  Ymax  STEP  ABS<Ytic> 

1468  MOVE  Xorg-Xfudge, I 
1470  LABEL  USING  1360;  I 
1488  NEXT  I 
1490  SUBEXIT 


"FREDIS"  (Cont.) 
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10 

DIM  FC1024) 

!  “PEAKS" 

20 

F$= “ T 1 1 V22 : F " 

30 

L=  1 9 

40 

PRINTER  IS  0 

50 

PRINT  F$, “PERKS  ABOVE '* , L 

60 

PRINT 

70 

PRINT  "ROW" , “COLUMN" , “POSITION", 

“AMPLITUDE" 

SO 

PRINTER  IS  16 

90 

ASSIGN  #1  TO  F $ 

100 

ON  END  #1  GOTO  120 

1  10 

READ  #1 ;F<*> 

120 

FOR  M  = 1  TO  32 

!  OUTPUT  DATA 

1  3  0 

C  =  0 

140 

FOR  N  =  M  TO  1023+M  STEP  32 

150 

C-C  +  l 

160 

IF  F < N ) > L  THEN  200 

170 

IMAGE  «, 4D, "  “ 

180 

PRINT  USING  170;FCN> 

190 

GOTO  250 

200 

IMAGE  # ,  4  D ,  "  *  " 

210 

PRINT  USING  200; F<N> 

220 

PRINTER  IS  0 

230 

PRINT  M,C,N,F<N> 

240 

PRINTER  IS  16 

250 

NEXT  N 

260 

PRINT  LINO) 

270 

NEXT  M 

!  NEXT  M 

280 

PRINTER  IS  0 

290 

PRINT  LINO) 

300 

END 

"PEAKS'* 
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IQ  DIM  F < 1 024  >  !  "DAT AST" 

20  INPUT  "AMPL I  TUBES  HIGHER  THEN" ,L 
30  FOR  J=10  TO  58 

40  PRINT  LINO) 

50  PRINT  “FRAME  NUMBER %J 

60  PRINT  LINO) 

70  S*=VAL *<J> 

80  ASSIGN  #1  TO  ,,TllV,,$cSn:,,:FM 

90  ON  END  #1  GOTO  110 

100  READ  #  1 ;  F  <  * ) 

110  FOR  M  = 1  TO  32  !  OUTPUT  DATA  <M> 

120  FOR  N=M  TO  1023+M  STEP  32 

130  IF  F  (  N  >  >  L  THEN  170 

140  IMAGE  #,4D, ”  M 

150  PRINT  USING  14Q;F(H> 

160  GOTO  190 

170  IMAGE  #, 4D, M*“ 

180  PRINT  USING  170;F<N> 

190  NEXT  N 

208  PRINT  LINCl) 

210  NEXT  M  !  NEXT  M 

220  NEXT  J 

230  END 


"DATAST" 


30 


£1 


lA&AwtfaflMihit 
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lCi  INPUT  “LEVEL 11  ,  L 

20  FOR  1=10  TO  50 

30  DIM  F < 1024) 

4d  1$=VAL$<I> 

50  vi*"Tiiru$r  :F" 

60  ASSIGN  #1  Tu  V# 

70  ON  END  #1  GOTO  90 

80  READ  # 1 ; F < * > 

90  GRAPHICS 

100  LINE  TYPE  1 

110  SCALE  0,12,0,15 

120  FOR  M= 1 5  TO  25 

130  C=0 

140  FOR  N=M  TO  1023+M  STEP  32 
150  C=C+ 1 

160  IF  F ( N ) < L  THEN  200 

170  MOVE  M- 1 5 , C 

180  LABEL  F ( N ) 

190  IF  C= 1 3  THEN  210 

200  NEXT  N 

210  NEXT  M 

220  WAIT  3000 

230  GCLEAR 

240  SCALE  0, 13,0,230 

250  FOR  M= 1 9  TO  22  !  PLOT  COLUMNS 

260  WAIT  200 

270  C=0 

280  GCLEAR 

290  MOVE  0 , L  !  PLOT  LEVEL 

300  LINE  TYPE  4 

310  DRAW  13, L 


!  “FRAME “ 


!  DATA  FROM  DISC 


!  LABEL  AMPLITUDE 


"FRAME" 
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320  LIME  TYPE  1 

330  MOVE  1,150 

348  LABEL  "COLUMN" 

350  MOVE  0,0 

360  FOR  N=M  TO  1023+M  STEP  32 

370  C  =  C+ 1 

380  PLOT  C-1,F<N> 

390  PLOT  C,F<N) 

400  NEXT  N 

410  WAIT  500 

420  MOVE  0,0 

430  NEXT  M 

440  GCLEAR  ,  pL0T  rows 

450  MOVE  0,0 

480  POP  N 1 =239  TO  335  STEP  32  ,  pL0T  LEvEL 

470  MOVE  0,L 

480  LINE  TYPE  6 

490  DRAW  13, L 

500  LINE  TYPE  1 

510  MOVE  1, 150 

520  LABEL  “ROW" 

530  MOVE  0,0 
540  POP  (1*1  TO  13 
550  N=N 1 +M 
560  PLOT  M-1,F<N> 

570  PLOT  M,FtN> 

580  NEXT  M 

590  MOVE  0,0 

600  WAIT  500 

610  GCLEAR 
620  NEXT  N1 

630  NEXT  1 

64W  END 


"FRAME"  (Cont.) 
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